Sensitivity of double resonance alignment magnetometers 
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We present an experimental study of the intrinsic magnetometric sensitivity of an optical/rf- 
frequency double resonance magnetometer in which linearly polarized laser light is used in the optical 
pumping and detection processes. We show that a semi-empirical model of the magnetometer can be 
used to describe the magnetic resonance spectra. Then, we present an efficient method to predict the 
optimum operating point of the magnetometer, i.e., the light power and rf Rabi frequency providing 
maximum magnetometric sensitivity. Finally, we apply the method to investigate the evolution of the 
optimum operating point with temperature. The method is very efficient to determine relaxation 
rates and thus allowed us to determine the three coUisional disalignment cross sections for the 
components of the alignment tensor. Both first and second harmonic signals from the magnetometer 
are considered and compared. 
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I. INTRODUCTION 

Our group develops optically pumped alkali vapor 
magnetometers (0PM) for both applied [l| and funda- 
mental [3] research. The diverse requirements of these 
demanding applications, in terms of sensitivity, spatial 
resolution, scalability, and measurement bandwidth, war- 
rants investigation of new 0PM schemes. An interesting 
and promising avenue is the use of atomic alignment in- 
stead of orientation to probe the external magnetic field. 
We will refer to an 0PM based on atomic alignment 
as DRAM (double resonance alignment magnetometer), 
while we will speak of DROM (double resonance orien- 
tation magnetometer) when the magnetization has the 
symmetry of an atomic orientation. Recently, our group 
presented both theoretical 3] and experimental [4l] in- 
vestigations of the magnetic resonance spectra produced 
in a cesium vapor in which an alignment is created and 
detected by a single linearly-polarized laser beam. 

Of direct importance for us, the DRAM scheme has a 
more flexible geometry than the well-known DROM M^,- 
configuration [Sj]. For maximal sensitivity, the DROM 
scheme requires a 45 deg. angle between the laser beam 
and the magnetic field [5|, limiting applications calling 
for a compact arrangement of multiple sensors. In mul- 
tichannel devices, as required for cardiomagnetic mea- 
surements for example, the DRAM method offers the 
advantage that the laser beam can be oriented either par- 
allel or perpendicular to the offset field without loss of 
sensitivity. 

The line shapes of the second harmonic DRAM signal 
have significantly narrower linewidths than the DROM 
signal under identical conditions. Narrow linewidths sup- 
press systematic errors in optical magnetometers, visible 
as long term baseline drifts, and potentially increase the 
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magnetometric sensitivity. This means that a DRAM 
could lead to a higher magnetometric sensitivity than a 
DROM for equal signal to noise ratio. 

Moreover, in DROM devices the interaction of the 
atoms with the circularly polarized laser light leads to 
an M dependent energy shift of the Zeeman hyperfine 
components when the laser frequency is not centered on 
the optical resonance line, the so-called light shift 0, [ll- 
In that case, the effect of laser power and frequency 
changes is indistinguishable from the effect of magnetic 
field changes, thus limiting the magnetometric perfor- 
mance and introducing systematic uncertainties on the 
determination of the absolute value of the field. In the 
DRAM, the linearly polarized light produces a light shift 
depending on M^, which does not have the same char- 
acteristics as a magnetic field Zeeman interaction. The 
A/^ shift broadens the magnetic resonance line, thereby 
slightly reducing the magnetic sensitivity, but it will not 
shift the resonance frequency. The absence of linear M- 
dependent systematic resonance shifts make the DRAM 
an attractive magnetometer for precision experiments 
searching for M-dependent effects, such as electric dipole 
moment searches [9|. 

In this article, we present an experimental study of the 
magnetometric sensitivity of a double resonance align- 
ment magnetometer. The principle of the DRAM with 
its theoretical description is given in Sec. [TTl the experi- 
mental setup is described in Sec. IIIIl and the operational 
definition of the magnetometric sensitivity is introduced 
in Sec. lIVI Then, in Sec.|V]we show that a simple empir- 
ical extension of the DRAM model extends its validity to 
significantly higher laser powers. This extended model is 
used in Sec. lVIl where we develop a method to predict the 
optimum operating point of the DRAM based on physi- 
cal parameters extracted from specific measurements. Fi- 
nally, the method is applied in Sec. IVIII to determine the 
temperature dependence of the DRAM optimum operat- 
ing point, and the results obtained for different cells are 
compared in Sec. IVIIII 



2 




FIG. 1: Double resonance magnetometer geometry using lin- 
early polarized light. Here, fe is the direction of the linearly 
polarized laser beam. The rf field Bi (shown here at t = 0) 
rotates in a plane perpendicular to the static field Bo- The 
linear polarization vector e makes an angle 7 with the static 
field Bo, and the phase of Bi is characterized by a. 



II. DOUBLE RESONANCE ALIGNMENT 
MAGNETOMETER 

The geometry of a double resonance alignment magne- 
tometer is presented in Fig. [T] it is identical to the one 
described in Q. A linearly polarized laser beam, with 
polarization e inclined at angle 7 to the magnetic field 
to be measured, Bq, is used to create an atomic align- 
ment via optical pumping in a room temperature vapor 
of cesium atoms. This alignment precesses under the si- 
multaneous action of the static magnetic field Bq and a 
much weaker magnetic field Bi, called the rf field, rotat- 
ing at frequency lu in the plane perpendicular to Bq and 
driving the magnetic resonance transitions. Competition 
between relaxation, optical pumping, and magnetic res- 
onance produces a steady state in the rotating frame. 

The precession of the alignment generates modulations 
of the absorption coefficient which create signals at both 
the fundamental (w) and the second harmonic (2uj) of the 
applied rf frequency u. The magnetic resonance signals 
Sui{t) and S2uj{t) are obtained here by monitoring the 
transmitted light power with a photodiode. The details 
of the calculation of Sui{t) and S2ui{t) are given in 
therefore only the most relevant equations needed for the 
magnetometric analysis are reproduced here. 

The magnetic resonance signals can be written as 

Suj (t) = Ao K (7) [ D^{5) cos {ujt - a) 

—A^{5) sa\{ujt — a)] , (la) 
^2^(t) = A /'-2c. (7) [ -A2u.{5) cos {2ujt - 2a) 

-D2u,{5) sin{2ujt-2a)] , (lb) 

where Aq is the alignment, defined in Q, produced by 
the optical pumping. The angular dependencies of the 
first and second harmonic signals are given by 
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(2 sin 27 + 3 sin 47) 



(2a) 



where 7 is the angle between the light polarization and 
the static field Bq. The first and second harmonic sig- 
nals have both absorptive, A^{5), A2^{5), and dispersive, 
D^{S), D2u]{5), components in their line shapes, given by 
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with a resonance denominator, 
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(l-4cos27-f 3cos47), (2b) 



In Eqs. (m and ([4]), ui — jpBi is the Rabi frequency 
of the rf field where is the gyromagnetic ratio of the 
ground state hyperfine level F. The detuning 6 = uj — loq 
is the difference between the rf frequency tu and the Lar- 
mor frequency ujq = ^fBq, and Fq, Fi, F2 are alignment 
relaxation rates. More precisely, the DRAM model Q 
calculates the evolution of the alignment multipole mo- 
ments TO2,g via a density matrix approach (for a general 
introduction to the use of multipole moments in the den- 
sity matrix formalism, see (To|). The moments m2.q are 
defined with respect to a quantization axis aligned with 
Bo and relax with rates T\q\. In practice, both the ab- 
sorptive A ((5), A2uj[5), and dispersive, D^^{5), D2ui{5), 
components of the signals can be used to measure the 
magnetic field, and can be isolated by phase-sensitive de- 
tection of the transmission signals Sui (t) and S2uj (t) . 

The effect of 7, the angle between the linear polar- 
ization vector and the static magnetic field, is contained 
in the functions h^(^) and /i2w(7)- The first harmonic 
signal is maximized for 7 = 25.5 deg and the second har- 
monic signal for 7 = 90 deg. Thus we distinguish between 
two realizations of the DRAM: 

1. the first harmonic DRAM, choosing 7 — 25.5 deg 
and measuring Ai^{5), D^{5), and 

2. the second harmonic DRAM, choosing 7 = 90 deg 
and measuring A2ui{6), D2ui{5)- 

Since the line shapes are different, we expect the two 
realizations of the DRAM to result in distinct optimum 
operating points and magnetometric sensitivities. 



III. EXPERIMENTAL SETUP 

The experimental setup used for the optimization pro- 
cedure is shown in Fig. [51 A Pyrex cell, paraffin-coated 



3 



polarizer 



XJ2 



Linear 
polarization 







.Bo 














cell 




1 




coils 



PD 



i(t) 



□FB 
laser 



DAVLL 



Ramp 




rf 




generator 




generator 










to 


— ► to tine rf coils 


Personal 




Lock In 




Computer 


QoiS) 


Amplifier 


Sir) = Ri(t) 



FIG. 2: Experimental setup: A cell containing Cs vapor was 
mounted inside a 3-axis Helmholtz coil array all placed inside 
a three-layer mu-metal shield. The polarization angle 7, mea- 
sured with respect to the offset field Bo, was set by a linear 
polarizer followed by a half-wave plate (A/2) located outside 
the shield for ease of access. The laser light traversed the cell 
and was converted to a current by a nonmagnetic photodiode 
(PD). All details can be found in the text. 



for spin relaxation suppression and evacuated except for 
an atomic cesium vapor in thermal equilibrium with a 
metal droplet, provided the paramagnetic atom sample. 
The cell was isolated from ambient magnetic fields by a 
three-layer mu-metal shield (ID 300 mm, length 580 mm, 
OD — 590 mm). Inside the shield, a primary pair of 
Helmholtz coils produced a static magnetic field Bo of 
about 3 /iT perpendicular to the light propagation direc- 
tion. Additional orthogonal pairs of Helmholtz coils (only 
one pair is shown in Fig. [J) were used to suppress residual 
fields and gradients. An rf magnetic field Bi, rotating at 
approximately 10 kHz in the plane perpendicular to the 
static magnetic field, was created by a set of two pairs 
of Helmholtz coils, wound on the same supports as the 
static field coils. All internal structural components were 
made from nonmagnetic materials. 

The laser beam used to pump and probe the atomic 
vapor confined in the cell was generated by a distributed 
feedback (DFB) diode laser, with a wavelength of 894 nm, 
stabilized to the 6Si/2,Fg = 4 ^ 6Pi/2,i^e — 3 hyper- 
fine transition by means of a dichroic atomic vapor laser 
lock (DAVLL) . A linear polarizer followed by a half- 
wave plate prepared linearly polarized light of adjustable 
orientation 7 with respect to Bq. The residual circu- 
lar polarization contamination was measured to be less 
than 1%. A nonmagnetic photodiode, followed by a low- 
noise transimpedance amplifier, detected the light power 
transmitted through the cell. The resulting signal was 
analyzed by a lock-in amplifier tuned either to the first 
or second harmonic of the rf frequency, depending on the 
DRAM configuration under study (cf. Sec.|TT|. A com- 
puter recorded magnetic resonance spectra by initiating 



the rf frequency sweep and simultaneously recording the 
in-phase and quadrature signals from the lock-in ampli- 
fier. The computer also controlled the light and rf power 
delivered to the cell, and measured the total light power 
on the photodiode as well as the temperature of the 
apparatus. The system was thus automated and could 
make measurements of the magnetic resonance signals 
for ranges of light and rf powers. Forced air heating was 
used to make temperature changes to the system, changes 
that were slow with respect to the time needed to record 
one spectrum. 

In practice, lock-in detection of the signals given by 
Eqs. HI) with respect to the rf frequency uj adds a phase 
4>i (selectable in the lock-in amplifier) and a small pick-up 
signal p(i 2)(yi,D) (smaller than 1% of the signal at maxi- 
mum) to each of the line shapes given by Eqs. ^ \^ . Due 
to (pi, the in-phase and quadrature spectra are, in gen- 
eral, a mixture of dispersive and absorptive line shapes. 
Demodulation of the signal Eq. ((Ta|, yields expressions 
used to fit the recorded in-phase and quadrature spectra 



lU^) = gUPL)hu{l) [ {D^{S) + Pid) cos ia+(f>i) 

+ {A^{6)+pia) sin {a + (pi)], (5a) 

QU^) = 9u{Pl)K{i) [ {A^{5) + Via) cos {a + 4>i) 

~{D^{6) +pid) sin {a + (pi)] . (5b) 

The guiiPh) factor is used here to contain not only ampli- 
fier gain factors, but also the alignment ^0 and any light 
power, Pl, dependencies. A similar mix of A2uj{5) and 
D2LJ {5) was used for the second harmonic signal given by 
Eq. HEI. 

Typical measured signals for the first and second har- 
monic magnetic resonance spectra are presented in Fig. [31 
together with fits of the theoretical line shapes given 
by Eqs. ([S])-®. All four curves are fitted simultane- 
ously with one set of relaxation rates; for the presented 
data To = 27r X 1.64(2) Hz, Fi = 27r x 2.93(2) Hz and 
F2 = 27r X 3.08(2) Hz. Detailed information on the fitting 
procedure is found in The excellent quality of the 
fits allow us to extract the amplitude g, the relaxation 
rates F^, and the Rabi frequency wi from a single set of 
double resonance spectra. For that reason no calibration 
of the rf-coils is needed. This is an advantage compared 
to the DROM where the Rabi frequency and the longi- 
tudinal relaxation rate are correlated to the point where 
they cannot be individually extracted from measured line 
shapes. 



IV. MAGNETOMETRIC SENSITIVITY 

The dispersive magnetic resonance line shapes given 
by Eqs. pa|) and ([Sc]) have a linear dependence on the 
detuning 5 = lo — ^fBq at the center of the resonance. 
By proper choice of (pi, the quadrature signal, Eq. (j5b[) . 
can be made completely dispersive, giving direct access 
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FIG. 3: (Color online) Measurements (circles) of the in-phase 
and quadrature magnetic resonance signals detected as am- 
plitude modulations of the transmitted light power, (a) First 
harmonic signals, (b) Second harmonic signals. The statis- 
tical uncertainty on each data point is represented by the 
symbol size. The solid lines are fits of the theoretical line 
shapes given by Eqs. ©-(IS]). 



to the linear zero crossing of the resonance, 

QUBo) - 9UPL)hUl) - ifBo) 



(6) 



with a similar expression for the 2w resonance. At con- 
stant w, can be used as a magnetometer signal for a 
limited range of magnetic field strengths \u} — ^fBq\ <C F, 
where F is the resonance linewidth. In that range, a 
change of i?o by a small Ai?o can be measured as a change 
of by the amount 



dBn 



(7) 



BQ=u/if 



The slope t„ is obtained from fits of the dispersive ex- 
perimental magnetic resonance line shape (Fig. [3]) using 
the relation 



dQu 
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Again, similar relations were used for the 2a; signals. 

The noise equivalent magnetic field, or NEM, repre- 
sents the noise limit on the derived value of B given the 



noise in Q. The total noise in Q has contributions from 
external magnetic field fluctuations and from all sources 
of technical noise, like laser intensity and frequency noise 
(converted to intensity noise by the atomic vapor), elec- 
tronic noise, and so on. All technical noises can, in princi- 
ple, be reduced until the system reaches the fundamental 
limit arising from the photocurrent shot noise. There- 
fore, we use the shot noise limited NEM as the measure 
for compariiig the performance of different magnetomet- 
ric schemes 

The root spectral density of the photocurrent shot 
noise is given by 



Ps = R y/ 2elDC , 



(9) 



where R is the transimpedance gain of the current ampli- 
fier, e the electron charge, and Idc the DC photocurrent. 
Given ps, the NEM can be expressed as a root spectral 
density of field fluctuations by inverting Eq. ([7|) 



NEM== 



PS 

\U 



(10) 



Since ps was evaluated from a measurement of the pho- 
tocurrent before the lock-in amplifier, the internal gain 
correction of the lock-in was used to give a measurement 
of t„ usable in Eq. (fTO|) . 

The goal of this study was to find the optimum DRAM 
operating parameters, and uji, yielding maximum 
magnetometric sensitivity, i.e., minimal NEM. 



V. EMPIRICAL EXTENSION OF THE DRAM 
MODEL 

As discussed in [1, the analytical expressions for 
the DRAM model [Eqs. dl-llll)] are valid for low laser 
power only, however, empirical formulas were presented 
modeling the light power dependence of the relaxation 
rates and of the global amplitude factors of the DRAM 
signals. Here, we present improved empirical formulas 
extending the DRAM model to even higher light powers, 
our goal being to cover the power domain that must be 
explored while optimizing the magnetometer. 

The following empirical formula successfully represents 
the laser power dependence of the first harmonic signal 



Pi 



{Psi+Pl){Ps2 + Pl) 



(11) 



where C is a constant and Psi, Ps2 are experimentally 
determined saturation powers for which we currently 
have no rigorous model in terms of fundamental phys- 
ical constants and processes. A similar formula applies 
to the second harmonic amplitude, but requires different 
values for both the constant and the saturation powers. 
The model reflects the expectation that both the creation 
of alignment as well as the ability to probe the alignment 
will saturate with increasing power. 



5 



In a similar way, the Pl dependence of the relaxation 
rates has been modeled by a power series and good agree- 
ment with the measured data was found using a second 
order polynomial for each rate 



ri(PL) 



Too 
Tio 



L J 



L 7 



T20 + a2PL + hPl 



(12a) 
(12b) 
(12c) 



We call the following parameter set the extended DRAM 
model parameters, 

{C,Psi,Ps2,'^oo,oiQ,l3Q,TiQ,ai,l3i,T2Q,a2,l32} (13) 

and note that they have to be determined experimentally. 
For that purpose, we have measured a series of double res- 
onance spectra as a function of laser power, and extracted 
the amplitude and relaxation rates from the simultane- 
ous fits, using common parameters, of the theoretical line 
shapes given by Eqs. Q-jS]) to the experimental data, 
as explained in [J. The measurements were made sep- 
arately for the first harmonic, with 7 = 25.5 deg, and 
for the second harmonic, with 7 = 90 deg. The results 
are presented in Fig. UHS]) for the first (second) harmonic 
signals. 

For determining C, -Pgi, and Ps2, the empirical model 
of Eq. ([TT|) was fitted to the amplitude data, and the re- 
sulting fits are displayed as solid lines in the upper graphs 
of Figs, m and [5] To find the remaining extended pa- 
rameters. Too, Fio, r2o, ao, ai, a2, /3o, Pi, and /32, the 
empirical model of Eqs. was fitted to the Pl depen- 
dence of the measured relaxation rates, and the resulting 
fits are displayed as solid lines in the lower graphs of 
Figs. [3] and [5] Clearly, the extended model accurately 
represents the data over the whole range of light powers 
investigated. 

Table U summarizes the extended DRAM model pa- 
rameters for both the first and the second harmonic sig- 
nals. The expectation, based on the cylindrical symmetry 
of the physical system, that Fio should equal is not 
refiected in the data, but the discussion of this will be 
delayed until Sec. En] 



VI. PREDICTION OF THE DRAM OPTIMUM 
OPERATING POINT 

A. Description of the method 

The optimum operating point of a DRAM can be pre- 
dicted from the measured extended DRAM parameters 
presented in the previous section. Here, the optimum 
operating point refers to the laser power P^ and Rabi 
frequency loi which minimize the intrinsic NEM defined 
in Eq. pUj) . In that equation, the photocurrent shot noise 
Ps is calculated from the DC photocurrent using Eq. ([9]), 
and the on-resonance slope of the magnetometer signal 
is calculated from the derivative of the dispersive compo- 
nent of the resonance spectra, see Eqs. ^ and ([5]). By 
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FIG. 4: (a) First harmonic DRAM signal amplitude as a func- 
tion of laser power, (b) Relaxation rates as a function of laser 
power. Points are measured values, extracted from the fit 
of the DRAM model [Eqs. ((3)l] to the experimental magnetic 
resonance spectra. The solid lines are fits of the extended 
DRAM model [Eqs. (fTTjl and (fT2|) ] to the experimental data. 
The data were measured at room temperature from a first 
harmonic DRAM with 7 = 25.5 deg, oji = 27r x 8.3 Hz. 



direct differentiation, we obtain 



(FiFs + ujD [FoFiFs -f (Fo + 3F2) u^l] 
for the slope of the first harmonic signal, and 



d5 



dQ 



2uj 



dS 



(5=0 



g2^(ii)fe2c.(7)ro(2Fi+F2)c^f 
(F1F2 + Lol) [F0F1F2 + (Fo + 3F2) Lol] 



(14) 



(15) 



for the slope of the second harmonic signal. Combin- 
ing the above with the power scaling model of Eqs. (fTT|) 
and (fT2)) and using the result in Eq. (fTO)) . the intrinsic 
NEM as a function of laser power P^ and Rabi frequency 
uji is found. 

This NEM function has been calculated for the ex- 
tended DRAM model parameters given in Table [J The 
resulting contour plots of NEM as a function of Pi and 
ui are presented in the upper graph of Fig. [5] for the first 
harmonic DRAM, and in the upper graph of Fig. [7] for 
the second harmonic DRAM. Both graphs show a clear 
optimum point where the NEM is minimum. Table |TT] 
lists the coordinates of these optimum points, together 
with the corresponding NEM value. 
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FIG. 5: (a) Second harmonic DRAM signal amplitude as a 
function of laser power, (b) Relaxation rates as a function of 
laser power. Points are measured values extracted from the fit 
of the DRAM model [Eqs. ((3)l] to the experimental magnetic 
resonance spectra. The solid lines are fits of the extended 
DRAM model [Eqs. (fTTj) and ([T2)) ] to the experimental data. 
The data were measured at room temperature from a second 
harmonic DRAM with 7 = 90deg, ui = 2n x 8.3 Hz. 
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TABLE I: Experimental values of the extended DRAM model 
parameters extracted from the fits to the experimental data 
presented in Figs. |4] and [S] Only statistical uncertainties are 
shown. See text for details. 



Fit pcLrctmeters 


■-y — 25.5 dsg 




c 


3.1(3) V 


1.0(2) V 


Psi 


1.7(2) fiW 


1.2(2) fiW 


PS2 


24(3) /iW 


11(3) fiW 




1.86(1) Hz 


2.03(3) Hz 


rio/27r 


3.27(1) Hz 


3.34(2) Hz 


r2o/27r 


3.58(5) Hz 


3.31(1) Hz 


Qo/27r 


0.492(4) Hz//iW 


1.21(1) Rz/fiW 


ai/27r 


0.934(7) Hz//^W 


1.33(1) Hz/mW 


Q2/27r 


1.55(3) Rz/fiW 


0.946(4) Hz/aiW 


/3o/27r 


-5.6(5) mRz/nW^ 


-28(2) mHz//iW2 


Pl/27V 


-8.9(9) mHz//iW2 


-24(1) mHz/z^W^ 


02/271 


-26(3) mRz/fiW^ 


-13.4(5) mUz/fiW^ 



TABLE II: Theory predictions of laser power _FL and rf field 
Rabi frequency wi minimizing the NEM compared to the ex- 
perimental best values. The calculations used the empirical 
extension of the DRAM model from Sec. |V] 



DRAM 


Optimum Ft 


Optimum uji/2n 


NEM 


scheme 


(mW) 


(Hz) 


(fT/x/H^) 


Theo. Ilu 


5.4 


2.34 


35.5 


Expt. lu! 


5.1(2) 


2.40(5) 


35.7(7) 


Theo. 2uJ 


4.2 


5.4 


32.8 


Expt. 2lu 


4.5(2) 


5.6(1) 


32.6(6) 



For the first harmonic DRAM, Fig. [6] shows a diagonal 
valley along oji = r2(ii)/V2 where the NEM IS maxi- 
mized (i.e., poor sensitivity). There, the NEM goes to 
infinity due to the onset of the narrow spectral feature 
(discussed in detail in [^) appearing on the dispersive 
component of the resonance spectra, cf. Fig.[3la), reduc- 
ing the slope to zero. 



B. Verification of the method 

The apparatus described in Sec. IIIII was used to mea- 
sure the on-resonance slope of the dispersive magnetic 
resonance signal. Then, that slope was inserted in 
Eq. (|10p to determine the experimental value of the in- 
trinsic NEM. We repeated the measurement on a 18 x 14 
grid of uJi and Pl values for the first harmonic DRAM, 
and a 9 X 15 value grid for the second harmonic DRAM. 
These measured points are shown as dots in the upper 
graphs of Figs. inland [71 For all measured points, the dif- 
ference between the NEM predicted from the extended 
model and the measured value has been determined, and 



the distribution of the relative difference is plotted in 
the lower graph of Fig. [6l for the first harmonic DRAM, 
and in the lower graph of Fig. [71 for the second harmonic 
DRAM. Within the experimental uncertainty, there are 
no significant differences between the measurements and 
the predictions. 

The experimental optimum operating points, where 
the measured NEM is minimized, was found, and the re- 
sults are shown in Table HIl Note that the optimal laser 
power is nearly the same for both first and second har- 
monic DRAMs. The Rabi frequency required to optimize 
the 2uj NEM is more than twice that of the Ilu NEM. The 
minimum NEM is slightly lower for the second harmonic 
signal. Table [Hi compares the measured values with the 
predicted values calculated using the extended DRAM 
model. The agreement is very good, in particular for the 
NEM values. This means that given the relaxation rates 
and saturation powers, the optimum point can be pre- 
dicted with precision of 5% using the extended DRAM 
model. 
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FIG. 6: First harmonic DRAM, 7 = 25.5 deg. The upper 
graph is a contour plot of the NEM as a function of Pl and 
uji/2tt. The NEM values were calculated using the method 
developed in Sec. I VI A1 The cross indicates the position where 
the NEM is minimum (see Table lll)l . The contour lines start 
at 40 fT/^Hz and are spaced by 10 fT/VlIz. The dots indi- 
cate the points in parameter space where the NEM has been 
measured, cf. Sec. I VI Bl The distribution of the relative dif- 
ference between calculated and measured values is shown in 
the lower graph. 
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FIG. 7: Second harmonic DRAM, 7 = 90 deg. The upper 
graph is a contour plot of the NEM as a function of Ft and 
uji/2tt. The NEM values were calculated using the method 
developed in Sec. IWAI The cross indicates the position where 
the NEM is minimum (see Table |TI|. The contour lines start 
at 40 f^/^/H^ and are spaced by 10 FT/VhS. The dots indi- 
cate the points in parameter space where the NEM has been 
measured, cf. Sec. IVIBI The distribution of the relative dif- 
ference between calculated and measured values is shown in 
the lower graph. 



C. Advantage of the method 

The automated experimental determination of the op- 
timum operating point of a DRAM, for a given temper- 
ature, can take several tens of hours. Indeed, that was 
the case for the NEM measurements over the grid of Pl 
and uJi values presented above. 

By contrast, to find the optimum operating point based 
on the prediction of the extended DRAM model parame- 
ters requires only measurements as a function of Pl , since 
the uji dependence of the magnetic resonance spectra is 
perfectly described by the DRAM model presented in [^] . 
Thus, the measurement time needed for finding the op- 
timum can be reduced by one order of magnitude when 
using the above method to predict the optimum operat- 
ing point instead of exploring the whole bidimensional 
parameter space. This is particularly useful to make a 



rapid characterizing of the quahty of coated Cs cells, and 
when studying the optimum point as a function of tem- 
perature, the topic of the next section. 



VII. DRAM OPTIMUM POINT EVOLUTION 
WITH TEMPERATURE 

The atomic vapor density and atom velocity distribu- 
tion (and hence the interatomic and waU coUision rates), 
as well as the relaxation probability during individual 
wall collisions, depend on temperature. The alignment 
relaxation rates. To, Fi, and T2, depend in a nontriv- 
ial way on aU those parameters because of three main 
contributing processes, namely coUisional spin-exchange, 
wall-collision electron-spin randomization, and the reser- 
voir effect [1^ [1^ . Temperature thus has an impor- 
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FIG. 8: Relaxation rates as a function of laser power, mea- 
sured at two different temperatures. The empty symbols 
(□,A,0) represent ro,ri,r2 measured at T = 25°C. The filled 
symbols (BjAj^) represent Io,ri,l2 measured at T = 38°C. 
The solid and dashed lines are a fit of the extended DRAM 
model [Eqs. p2|l ] to the experimental data, (a) First har- 
monic DRAM with 7 = 25.5 deg. (b) Second harmonic 
DRAM with 7 = 90 deg. 



tant influence on the magnetometric sensitivity. Unfor- 
tunately, the influence is hard to model and so is worth 
measuring. 

The method developed in Sec.|Vl] — predicting the op- 
timum point from a measurement of the extended DRAM 
model parameters (cf. Sec. |V]) — was used to investigate 
the temperature dependence of the optimum DRAM op- 
erating point. Multiple measurements of the DRAM pa- 
rameters were made for temperatures between 20° C and 
40° C. As an illustration, the evolution of relaxation rates 
with temperature is shown in Fig. [8] where two measure- 
ments, at 25°C and 38°C, are presented. Even though 
the temperature evolution of the relaxation rates is non 
trivial, the quadratic model of Eqs. fits well to the 
experimental data and gives access to the relaxation rates 
Too J Tio, and at zero light power. These parameters 
have been measured and their temperature behavior is 
presented in Fig. [S] They all increase with tempera- 
ture and this is mainly related to the increase of the 
atomic vapor density. During setup, we observed that 
the difference between Fio and F20 can be decreased by 
improving the magnetic field homogeneity, and further 



FIG. 9: Temperature evolution of the relaxation rates extrap- 
olated to zero laser power. The symbols (l,^,^) represent 
loo, Fio, r2o. The solid lines are fits of Eq. (I16|l to the exper- 
imental data, (a) First harmonic DRAM with 7 = 25.5 deg. 
(b) Second harmonic DRAM with 7 = 90 deg. 

tests confirmed that the difference between Fio a-nd F20 in- 
creases with the square of the magnetic field inhomogene- 
ity Moreover, the residual difference observed in 

our experiment is compatible with an estimation of the 
residual magnetic field inhomogeneity. This leads us to 
conclude that, in principle, the two transverse relaxation 
rates Fio and F20 should be equal in a perfectly homoge- 
neous magnetic field. In Fig. [51 the solid lines are fits to 
the experimental data of the relaxation model given by 

Tio = n CTi Urci + Av^ e^''/^^ + Bv,^ + C v^^ . (16) 

On the right-hand side, the first term is the contribu- 
tion due to coUisional spin-exchange, it is proportional 
to the vapor density n{T), to the coUisional disalignment 
cross-section Ci, and to the atoms' mean relative velocity 
v,-ai{T) = y/WkT/inM) where M is the "^Cs mass. The 
second term is the contribution due to wall-collisions: it 
is proportional to the rate of wall collisions, h ence to the 
atoms' mean velocity v^iT) — yJSkT / {-kM), and to the 
wall sticking time Ts ~ Toe^"/'^^ where tq ~ 10~^^ s, Ea 
is the adsorption energy and k is the Boltzmann constant. 
The third term is the contribution from the reservoir ef- 
fect, and it is proportional to the rate of wall collisions 
and therefore to I'm- Finally, the last term is the con- 
tribution due to magnetic field inhomogeneities, which is 
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TABLE III: CoUisional disalignment cross-sections Ui ob- 
tained from the fit of Eq. (|16p to tlie experimental data pre- 
sented in Fig. [9l The last column gives the average of the 
values obtained for the luj and 2ijJ DRAMs. 



Parameter 


Ilo dram 


2uj DRAM 


Average 


(Jo (cm^) 
(Ti (cm^) 
CT2 (cm^) 


0.5(6) X 10^" 
1.5(3) X 10-" 
1.5(3) X 10-" 


1.2(6) X 10"" 
2.1(9) X 10"" 
1.7(5) X 10"" 


0.9(4) X 10"" 
1.6(3) X 10-" 
1.6(3) X 10"" 



proportional to due to motional narrowing p^ . 

The vapor density n{T) is calculated from the cesium 
vapor pressure given in it is highly temperature 

dependent. As a consequence, over the range of tem- 
peratures investigated in this work, the collisional spin- 
exchange term represents the main contribution to the 
temperature behavior of relaxation rates, and all other 
terms are approximately linear. Therefore, the fit is 
able to determine the Cs~Cs collisions disalignment cross- 
sections fJi , but cannot distinguish the contributions from 
the other terms with reasonable uncertainties. In princi- 
ple, this can be improved by increasing the temperature 
range of the measurements, and would lead to a powerful 
method for the investigation of relaxation mechanisms. 
However, at present, experimental setup cannot reach 
the necessary temperatures and so, since it is beyond the 
scope of the present paper, such investigations will be the 
subject of future work. The values of ui extracted from 
the fits are summarized in Table IIIII 

The extended DRAM model parameter measurements 
were used to calculate the evolution of the optimum op- 
erating point of the DRAM (cf. Sec. NTK^ . The results 
are presented, as a function of temperature, in Fig.fTOlfor 
both the first and second harmonic DRAMs. A quadratic 
polynomial was fitted to the NEM data in order to de- 
termine the temperature of minimum NEM. The exper- 
imental parameters characterizing these optimum points 
are summarized in the first column of Table IIVI where 
we observe that the second harmonic DRAM is slightly 
more sensitive than the first harmonic DRAM. 



VIII. DISCUSSION 

The majority of the work presented herein has been 
performed using a single evacuated Cs cell (Cell 1 in Ta- 
ble [IV])- We applied the temperature NEM optimization 
procedure described in Sec. lVIll to two additional paraffin 
coated cesium cells and the resulting optimum parame- 
ters are presented in Table IIVI The intrinsic NEM for 
the second harmonic DRAM is always smaller than for 
the first harmonic DRAM and the lowest NEM value of 
27.4 fT/-\/Hz was obtained using the cell produced by our 
group (Cell 1). We observe that the NEM values scale 
with the inverse of the volume of the cell and not with 
the volume to surface ratio. However, since the three 
cells do not have the same coating, this relation could be 
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FIG. 10; Evolution of the optimum NEM as a function of tem- 
perature, (a) First harmonic DRAM with 7 = 25.5 deg. (b) 
Second harmonic DRAM with 7 = 90 deg. These graphs were 
calculated from the measurement of the extended DRAM 
model parameters. 



accidental. 

Under identical conditions, the line shapes of the sec- 
ond harmonic DRAM signal are narrower than those of 
the DROM which, in principle, should lead to an im- 
provement of the sensitivity [3]. Previous work by our 
group found an intrinsic sensitivity of 10 fT/VHz for an 
optimized DROM using a 70 mm diameter Cs vapor cell 
in the so called M2;-configuration 18]. However, the ex- 
perimental setup used in the past was very different (dif- 
ferent cell size, offset field homogeneity, and magnetic 
shielding) and therefore a detailed comparative study is 
needed before drawing firm conclusions. 



IX. CONCLUSION 

In conclusion, we have presented an experimental 
study of the intrinsic magnetometric sensitivity of the 
double-resonance alignment magnetometer, showing that 
an empirical extension of the DRAM model can be used 
to describe the magnetic resonance spectra over a range 
of experimental parameters sufficient for optimizing the 
magnetometer. A model has been developed to predict 
the optimum operating point of the magnetometer, i.e., 
the value of experimental parameters for which the mag- 
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TABLE IV: Results of the temperature NEM optimization 
procedure (as described in Sec. IVII|I applied to three differ- 
ent paraffin coated cells. Cell 1 was produced by our group, 
it is spherical with a 28 mm diameter. Cells 2 and 3 were 
purchased from a Russian company. 





Cell 1 


Cell 2 


Cell 3 


Shape 


Spherical 


Cubic 


Cylindrical 


Volume V 


11.5 cm^ 


8.0 cm^* 


4.6 cm^ 


Surface S 


24.6 cm^ 


24.0 cm^ 


15.3 cm^ 


V/S 


0.467 cm 


0.333 cm 


0.301 cm 


First harmonic DRAM 


Optimum T 


31.2 °C 


35.8 °C 


35.5 °C 


Optimum Ft 


6.2 


12.1 fiW 


5.2 ^iW 


Optimum uii/2n 


2.9 Hz 


5.5 Hz 


4.8 Hz 


Intrinsic NEM 


28.6 fr/ySS 


45.7 FT/V^ 


74.7 FT/VhS 


Second harmonic DRAM 


Optimum T 


33.7 °C 


36.2 °C 


37.6 °C 


Optimum Ft 


6.0 


8.9 fiW 


4.7 


Optimum (jji/27r 


10.5 Hz 


14.1 Hz 


15.0 Hz 


Intrinsic NEM 


27.4 fr/VHi 


39.3 fr/%/Hi 


62.2 fl/VTh. 



netometric sensitivity is maximum. The method was ver- 
ified by comparing its results to a direct measurement of 
the optimum operating point. In contrast to the time 



consuming direct optimization involving many hours of 
testing in a two parameter space, our method decreases 
the time required to find the optimum operation point 
to half an hour. Finally, we used this method to investi- 
gate the evolution of the optimum operating point of the 
DRAM with temperature, showing that the magnetomet- 
ric sensitivity reaches an optimum of 27.4 fT/VHz for a 
temperature of 33. 7° C. Both the first harmonic and the 
second harmonic realizations of the magnetometer were 
explored and compared. The temperature dependence 
of the relaxation rates yielded measurements of the Cs- 
Cs collisional disalignment cross sections of the tensor 
alignment, and the method promises to be useful in the 
continued study of the relaxation processes over broader 
temperature ranges. 
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